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1 Introduction 

String theory (ST) and Vasiliev's equations (VE) [1, 2] are the only known examples of 
consistent theories of interacting higher-spin (HS) particles. 1 Although their current for- 
mulations provide mathematically elegant descriptions involving infinitely many auxiliary 
fields, some important aspects, as the number of derivatives involved in the cubic vertices 
or the possible (non-)local nature of higher-order interactions, are hidden. Indeed, the 
cubic vertices involving fields in the first Regge trajectory of the open bosonic string have 
been obtained only recently in [5, 6]. Leaving aside Chan-Patton factors, their illuminating 
form is given by 



where i,j = 1,2,3 label the Fock spaces associated to the interacting particles. From the 
latter expression, the s± — S2 — S3 interactions, as well as the corresponding coupling con- 
stants, can be easily extracted via a Taylor expansion of the exponential function. Besides 
reflecting the world-sheet (Gaussian integral) origin of ST, the latter fulfills all requirements 
dictated by the compatibility with both the string spectrum and the corresponding global 
symmetries. Hence, many key properties can be deduced by investigating the consistent 
cubic vertices. Concerning VE, one might ask what is the form of the cubic vertices and 
what can one learn from them. Given the analogy to ST, we expect to understand how 
the global HS symmetry constrains the massless s\ — S2 — S3 interactions making the entire 
spectrum of VE as a single immense multiplet. This question has been partly addressed in 
[7-9], where it has been shown that, starting from VE, the extraction of the cubic vertices 
requires infinitely many field redefinitions making the analysis very involved. 

On the other hand, moving from a top-down to a bottom-up viewpoint, one may ask 
oneself which HS cubic interactions lead to fully nonlinear HS theories and whether ST 
and VE are the only solutions or not. This question (called Gupta or Fronsdal program) 
can be tackled solving the Noether procedure for HS fields. The latter is a perturbative 
scheme (order by order in the number of fields) whose aim is to classify all consistent 
interactions starting from a given free theory. The first step of such procedure is to find 
out the most general couplings of three massive or massless HS particles in an arbitrary 
constant-curvature background. For the case of symmetric HS fields, this problem has 
been addressed in [10, 11], where, making use of the ambient-space formalism, all possible 
S1 — S2 — S3 interactions were provided. 2 

Notice that the aforementioned program is closely related to the classification of all 
consistent CFTs in arbitrary dimensions. More precisely, in the context of AdS/CFT, the 
Noether procedure seems to share many analogies with the problem of classifying all pos- 
sible consistent OPEs of HS operators Oi(x) . In turn, this is tantamount to enumerating 

1 See [3, 4] for some recent reviews. 

2 See [12, 13] for a frame-like approach to the problem of massless interactions, and [9, 14-19] for other 
works on HS cubic interactions in (A)dS. 
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all possible tensor structures of three-point functions (Oi(xi) Oj{x<i) O^x^)), from which 
all conformal blocks can be computed. 3 In this respect, it would be interesting to under- 
stand the dictionary between the Noether procedure requirements on the bulk side and 
the conformal symmetry on the boundary side. This perspective can possibly clarify the 
role of Lagrangian locality, usually assumed in the bulk, or of possible alternatives, and 
may provide a new look into the AdS/CFT correspondence itself. Moreover, it is worth 
mentioning that the CFT results (3-pt functions) can be obtained a priori from the AdS 
results (cubic vertices) attaching the Boundary-to-Bulk propagators to the vertices. 

In the present paper we present the construction of HS cubic interactions in (A)dS 
along the lines of [10, 11]. We shall show that this problem is equivalent to finding poly- 
nomial solutions of a rather simple set of linear PDEs. Each solution is in one-to-one 
correspondence to a consistent cubic interaction. Let us stress that, since the solution 
space is linear, an arbitrary linear combination of these cubic vertices is also consistent, 
leaving their relative coupling constants unfixed (at this order). The latter are constrained 
within the Noether procedure either by compatibility with the global symmetries of the free 
theory and/or by consistency of quartic interactions. Let us stress once again the connec- 
tion to the conformal bootstrap program which may entail key (still unclear) requirements 
dictated by the Noether procedure. 

The organization of the paper is the following: in Section 2 we introduce the Noether 
procedure which represents the main tool of our construction. In Section 3 we briefly 
review how to apply such a scheme to derive cubic interactions in flat space. Section 4 is 
devoted to the formulation of the free theories in the ambient space. The ambient-space 
action at the cubic level is discussed in Section 5. Finally, in Sections 6 and 7 we present 
the solution to the cubic order of the Noether procedure in (A)dS. 

2 Noether procedure 

The aim of the Noether procedure is to find all consistent (at least classically) inter- 
acting structures associated to a given set of particles, order by order in the number of 
fields. 4 In the case of massless spin 1 or spin 2 particles (An or hn V ), this corresponds to 
identifying the consistent interactions starting from Maxwell or Fierz-Pauli Lagrangians. 
Arbitrary vertices involving or hn U would mostly cause a propagation of unphysical 
DOFs, which, at the free level, are removed by the linear gauge symmetries: 5 m A^ = d^a 
and S m hfj, u = d^ u + d u ^ . Hence, a key condition for the consistency of the interacting 
theories is the existence of gauge symmetries which are nonlinear deformations of the linear 
ones. 

Let us consider an arbitrary set of gauge fields (p a (where a labels different fields) 
with free action S" (2) and linear gauge symmetries 5^(p a . The problem is to find the 
corresponding nonlinear action S together with the non-linear gauge symmetries 8ip a . For 

3 See [20-25] for the correlation functions of three conserved currents. 

4 See [26] for a detailed discussion. 
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this purpose, one can consider the following perturbative expansions: 



S = S (2) + S (4) + ■■■ , (2.1) 

5 <p a = 5 (( V + <5 (1 V a + 5 (2 V a + • • • , (2.2) 

where the superscript (n) denotes the number of fields involved. Taking the variation of 
the action (2.1) under the gauge transformation (2.2), one ends up with a system of gauge 
invariance conditions: 

5 W S {2) = 0, (2.3) 

5(0)5(3) +5 d) 5 (2) =Q) (2.4) 

5(0)5(4) + £(1)5(3) + 6 (2) S (2) =0] (2.5) 



The first equation implies the linear gauge invariance of the free theory. The second equa- 
tion is a condition for both the cubic interactions 5 (3) and the first-order gauge deformations 
5^tp a , and so on. Since the second term in (2.4) is proportional to the free EOM, condition 
(2.4) implies 

<5 (0) S (3) «0, (2.6) 

where ~ means equivalence up to the free EOM. Solving this equation one can identify all 
cubic interactions consistent with the linear gauge symmetries. In the case of massless spin 
1, one finds two independent interaction terms which schematically read 

5 (3) = A^ 1 - 1 JAAF + X 1 ^- 1 JFFF. (2.7) 

The first is the one-derivative YM vertex while the second is the three-derivative Born- 
Infeld one. In the massless spin 2 case, there are three independent interactions: 

5 (3) = A 2-2-2 J h Q hdh + X 2-2-2 J RR + ^-2-2 J RRRj (3.8) 

where the first is the two-derivative gravitational minimal coupling while the other two 
come from the expansions of (Riemann) 2 and (Riemann) 3 and involve four and six deriva- 
tives respectively. As one can see from these lower-spin examples, the general solutions 
to eq. (2.6) are s± — S2 — S3 vertices with different number of derivatives associated with 
the coupling constants An 1_S2_33 . It is worth noticing that these coupling constants are 
independent at this level. 

The next step consists in solving (2.4) for the first-order gauge transformations S w (p a 
associated to each solution 5* (3) found at the previous step. In the lower-spin cases, only 
the first vertices in (2.7) and (2.8) lead to nontrivial deformations: 

S W A = X 1 ^ 1 ' 1 Aa, 5 {1) h = X 2 2 2 - 2 (h 8£ - dh £) , (2.9) 

of the linear gauge transformations. The latter correspond to the standard non-Abelian 
YM gauge transformations and to diffeomorphisms respectively. Although not deforming 
the gauge symmetries, the remaining vertices can be completed to the full non-linear order 
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keeping consistency with the gauge transformations (2.9). They form the first elements of a 
class of higher-derivative gauge or diffeomorphism invariants, where the remaining elements 
appear at higher orders S {n) . General (HS) gauge theories present as well two types of 
cubic vertices: the ones deforming the linear gauge symmetries, and the ones giving rise to 
possible higher-derivative gauge invariants. Although the former define the deformed non- 
Abelian gauge algebra, the second are also relevant since they provide possible (quantum) 
counter-terms. Hence, if no independent non-deforming vertices survive at higher-orders, 
then no counter-terms would be available and the theory would be UV finite. This issue 
can be initially addressed solving eq. (2.5) which involves quartic interactions as well as 
the spectrum of the theory [27, 28]. 

So far we have only considered the case of gauge theories. Constructing interactions of 
massive HS fields also raises similar problems. Arbitrary interaction vertices would mostly 
violate the Fierz conditions resulting in the propagation of unphysical DOFs. One way 
to proceed is to introduce Stueckelberg symmetries into the massive theory and perform 
the Noether procedure. Another equivalent way, that is actually the one we use in the 
following, is to first construct the transverse and traceless parts of the vertices and then 
consistently complete them. 



3 Flat-space case 

In this section we consider flat-space cubic vertices since their construction reveals some 
of the key ideas also used in the (A)dS case. See the review [3] for an exhaustive list of 
references on the cubic interactions, and [29-32] for more recent developments. 

In order to deal with arbitrary HS fields, it is useful to introduce auxiliary variables 
which are the analogue of the string oscillators (a l ^ 1 o u^), and define the generating 
function: 

^(x,n) = i<..., s (x)u^---^. (3.1) 

Here, the superscript A labels different HS fields, and in the following we use A = a for 
massless fields and A = a for massive ones. In this notation, the most general form of the 
cubic vertices is 

S (3) = f d d xC AlA2A3 (d Xi ,d Ui )ip A ^x 1 ,u 1 )^(x 2 ,u 2 )'p A3 (x 3 ,u 3 ) Xi=x , (3.2) 

where i = 1,2,3. Different functions C AiA2 a 3 describe different vertices, embodying the 
coupling constants of the theory. Notice that C Ai a 2 a 3 plays the same role as the state 
| V3) in the BRST (String field theory) approach. Restricting the attention to the parity 
invariant interactions, the dependence of C Ai a 2 a 3 on the six vectors d x . and d Ui , is through 
the 21 (6+9+6) Lorentz scalars d x . ■ d Xj , d Ui ■ d Xj and d Ui ■ d Uj . For instance, a vertex of the 
form: 

^ vp xd"d v ^ a d T ^ x , (3.3) 

is encoded by 

C = (d Ul ■ d X2 ) 2 (d Ul ■ d U2 ) (d Ul ■ d U[1 ) (d U2 ■ d U2 ) (d U3 ■ d X3 ) . (3.4) 
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Notice that not all C Ai a 2 a 3 , s are physically distinguishable but there exist two kinds of 
ambiguities. The first is due to the triviality of total derivative terms (or integrations by 
parts). This ambiguity can be fixed by removing d Ui - d Xi _ 1 in terms of the other Lorentz 
scalar operators as 

d Ui ■ d Xi _ x = - d Ui ■ d Xi+1 - d Ui ■ d Xi + d Ui ■ d x [i ~ i + 3] . (3.5) 

The second ambiguity is related to the possibility of performing non-linear field redefinitions 
which can create fictive interaction terms. However, these vertices are all proportional to 
the linear EOM so that the corresponding ambiguity can be fixed by disregarding the on- 
shell vanishing vertices. This amounts to neglect the dependence of the function C Ai a 2 a 3 
on the d Xi - d x .'s, since the latter can be expressed as 

d Xi ■ d Xl _, = \ (d 2 Xi+ - d 2 Xi - dl^) +\d x - (d x . + - d Xi+1 ) , (3.6) 

and, up to EOM, the d% 's can be replaced by the d Ui - d Xi 's and the d 2, s. For instance, 
the (Fronsdal's) massless HS EOM reads 

[d 2 x - u ■ d x d u ■ d x + \ (u • d x f d 2 ] <p a (x, u) « o . (3.7) 

Taking into account the latter ambiguities, the vertex function C Ai a 2 a z can only depend 
on 12 (3+3+6) Lorentz scalars: d Ui ■ d Xi+1 , d u . ■ d Xi and d Ui ■ d Uj . It is worth noticing that 
at the cubic order, there are no non-local vertices since the d Xi - 9 x .'s have been removed 
while the other scalar operators can only enter with positive powers (otherwise the tensor 
contractions do not make sense). 

When a gauge field, say (p ai , enters the interaction, the function C aiA2 A- s is further 
constrained by the condition (2.6). In this notation, the linear HS gauge transformations, 
5 m vt-u = d(„ „ ^ , read 

5 m (p a (x,u) = u-d x e a (x,u). (3.8) 

Hence the cubic-order gauge consistency condition (2.6) gives 

[ C ai a 2 a 3 (d Ul ■ d Xi+1 , d Ui ■ d Xi , d Ui ■ d Uj ) , u x ■ d Xl ] wO, (3.9) 

where ~ means equivalence modulo the Fronsdal equation (3.7). In order to tackle the 
above equation, it is convenient to split the function C aiA2 A 3 into two parts: the one 
which does not involve any divergence, d Ui ■ d Xi , trace, d\. or auxiliary fields (we call it the 
transverse and traceless (TT) part), and the one which does (DTA part). Let us notice that 
the TT part is precisely what survives after eliminating unphysical DOF. Indeed, besides 
the mass-shell condition, the Fierz system: 

Fierz system : (Sg - m 2 ) ip A = , d u ■ d x y A = , d 2 u y A = , (3.10) 

involves also the transverse and traceless conditions. 5 Therefore, the TT part of the action 
plays a key role encoding the on-shell content of the theory. On the other hand, the 

5 When m = 0, one has to quotient the system by the gauge symmetries (3.8) with parameters e a 
satisfying the same conditions (3.10). 
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part containing divergences, traces or auxiliary fields vanishes after gauge fixing. 6 The 
next question is whether it is possible to determine the TT part of the vertex without 
using any information about the other part. From a physical point of view, this ought 
to be possible since the physical (on-shell) interactions cannot depend on the unphysical 
ones. 7 Concerning massive fields, the TT conditions already assure the propagation of the 
correct DOF and no further constraints has to be imposed on the TT parts of the cubic 
interactions. However, the remaining parts of the massive vertices are constrained in such 
a way that the Lagrangian EOM reduce to the Fierz system (3.10). 

In the following, we show how to determine the TT parts of C Ai a 2 a 3 from eq. (3.9). 
First, after removing all the ambiguities through eqs. (3.5,3.6), any functional F can be 
univocally written in terms of its TT part and the remaining part as F = [F] TT + [-F] D ta • 
Hence, eq. (2.6) can be split into two equations: 

[<W] TT «0, [<* (0) S (3) ] DTA «0, (3.11) 

where, henceforth, ~ means equivalence modulo the Fierz system (3.10). Second, as the 
gauge variations of divergences, traces or auxiliary fields are proportional to themselves 
up to c^.-terms: [5 (0) [5 (3) ] DTA ] TT ~ , the first equation in (3.11) gives an independent 
condition for the TT parts, [iS (3) ] TT , of the interactions: 

[w] TT = r { [^] TT + [^] DTA }] tt « r [^)] t j tt „ o. (3.i2) 

At this point, [5 <3) ] TT can be expressed as in eq. (3.2) through a function CJ^^^ (yi, z{) of 
6 variables: 

Hi = dui' d Xi+1 , z% = d Ui+1 - d u ._ 1 . (3.13) 

Then, assuming the first field to be massless, Ai = a\, eq. (3.12) gives a condition for 
^aiA 2 A 3 ana l°g° us to (3.9). Using the Leibniz rule, we obtain a rather simple differential 
equation: 

[m d Z3 - V3 d Z2 + \ {ml - ml) d yi ] = , (3.14) 

where rrii is the mass of the i-th field. When two or three massless fields are involved in 
the interactions, one has respectively two or three differential equations given by the cyclic 
permutations of eq. (3.14). 

Depending on the cases, the corresponding solutions CJ^^^ are constrained to depend 
on some of the y^s and the z^s only through the building blocks: 

9 = VlZl +V2Z2 +V3Z3 , (3.15) 
hi = y i+ i Vi-i + \ [ml - (m i+ i + mj_i) 2 ] m . (3.16) 

As an example, we consider the interactions of three massless HS fields where the consistent 
cubic interactions are encoded in an arbitrary function: 

C I?a 2 a 3 = £aia 2 a 3 0/1,2/2, 2/3, fiO • (3-17) 



6 Indeed the light-cone interaction vertices [33] can be obtained solely from the TT part by going to the 
light-cone gauge. 

7 The light-cone gauge approach is consistent in its own without calling for some additional conditions 
on its corresponding covariant off-shell description. 
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Leaving aside Chan-Paton factors, the latter can be expanded as 

min{si,S2,S3} 

ic= a™ 3 g n y s r n y s r n y s r n > (3-18) 

n=0 

where the A^ 1_S2_S3 's are independent coupling functions that ought to be fixed by the quest 
for consistency of higher order interactions. From the latter expression it is straightforward 
to see that the number of consistent couplings is min{si, S2, S3} + 1 , while the number of 
derivatives contained in each vertex is si + S2 + S3 — 2 n . In particular, focussing on the 
2 — 2 — 2 case, eq. (3.18) gives 

K = \r~ 2 9 2 + At 2 " 2 9 Vi V2 ys + Aq~ 2 ~ 2 y\ y 2 y 2 , (3.19) 
that exactly reproduces eq. (2.8). 

4 Ambient-space formalism for HS 

In order to address the HS interaction problem around an arbitrary constant-curvature 
background (i.e. (A)dS space), one can still rely on the Noether procedure introduced in 
Section 2. However, in this case the starting point are the HS free theories in (A)dS, where 
besides massive and massless particles, new types of particles (called partially-massless) 
appear [34]. 8 Moreover, the cubic interactions built on top of the free theories would 
involve (A)dS covariant derivatives whose non-commuting nature makes the construction 
cumbersome. Their commutators give rise to lower-derivative pieces proportional to the 
cosmological constant, making the vertices inhomogeneous in the number of derivatives. 
The ambient-space formalism proves to be a convenient tool in dealing with free (A)dS HS, 
and hence, it represents a natural framework in order to construct (A)dS cubic interactions. 
Furthermore, recently it has been intensively used in the context of Mellin amplitude in 
the computations of Witten diagrams [36, 37]. 

The ambient-space formalism [38, 39] consists in regarding the <i-dimensional (A)dS 
space as the hyper-surface X 2 = el? embedded into a (d + l)-dimensional flat-space. In 
our convention the ambient metric is r/ = (— , +, . . . , +), so that AdS (e = —1) is Euclidean 
while dS (e = 1) is Lorentzian. Focussing on the region eX 2 > 0, there exists an iso- 
morphism between symmetric tensor fields in (A)dS, v^...^ , and those in ambient space, 
^Afi—Afs ) satisfying the homogeneity and tangentiality (HT) conditions: 

Homogeneity : (X ■ d x - U ■ d n + 2 + n) $(X, U) = , 
Tangentiality : X ■ d v $(A, U) = . (4.1) 

Here we have used the auxiliary-variable notation for the ambient-space fields: 

$(X, U) = i $ Mr ..M s P0 U M ^ • • • U M * . (4.2) 

8 In the case of mixed-symmetry HS fields in AdS, even the notion of massless-ness changes with respect 
to the flat-space case [35]. 
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The degree of homogeneity \x parametrizes the (A)dS mass-squared term appearing in the 
(A)dS Lagrangian: 



m 



2 _ t 



L- 



[( l x-s + 2)( fJ ,-s-d + 3)-s], (4.3) 



so that = corresponds to the massless case. In the ambient-space formalism, the EOM 
of both massless and massive HS fields are given by the Fronsdal ones (3.7) after replacing 
(x, u) by (X, U) . Let us remind the reader that the concept of massless-ness in (A)dS 
is not related to the vanishing of the mass term but rather to the appearance of gauge 
symmetries. In fact, if one postulates the latter to be of the form: 

5 (0) <S>(X,U) = U-d x E(X,U), (4.4) 

then the compatibility with the HT conditions (4.1) alone restricts both the possible values 
of n and the normal(radial) components of E . In particular, when fi = 0, 1, . . . , s — 1, then 
there exist compatible higher-derivative gauge symmetries: 

5 m $(X, U) = {U- dxY +1 0(X, U) [E = (U-d x Tn], (4.5) 

with the gauge parameters O satisfying 

(X -d x -U -du - fi)n(X,U) = 0, X -dun(X,U) = 0. (4.6) 

On the other hand, for other values of /z, no gauge symmetries (in absence of auxiliary fields) 
are allowed, implying that the corresponding fields are massive. Notice that the massless 
field, fx = 0, is the first member of a class of representations where the other members, with 
[i = 1, 2, . . . , s — 1, are called partially-massless. However, partially-massless fields describe 
unitary representations only in dS. 

Before closing this section, let us discuss the flat limit from the ambient-space view- 
point. The latter consists first in translating the coordinate system as X M — > X M + L N M , 
where iV is a constant vector satisfying iV 2 = e , and second, in taking the L — > oo limit. 
As a result, the HT conditions (4.1) reduces to 

(N -d x -V^M)<Z>(X,U) = 0, N-du^(X,U) =0, (4.7) 

where the flat mass M is related to the (A)dS mass \i as 

v /z eM = - lim ^. (4.8) 

Notice that in this limit, all (A)dS representations become massless, while, in order to 
recover massive representations in flat space one should consider the /x — > oo limit. 



5 Ambient-space action 

In the previous section we have shown how to describe HS fields in (A)dS making use 
of the ambient-space language. In Section 6 we shall use the latter framework in order 
to solve the Noether procedure. For this purpose, one needs to know first of all how to 
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express the (A)dS action in terms of ambient-space quantities. As far as the Lagrangian 
is concerned, no subtleties arise since, together with the isomorphism between (A)dS and 
ambient-space fields, there is an analogous one between (A)dS-covariant derivatives V M 
and ambient-space ones 8 x m — (X M /X 2 ) X ■ dx- Hence, any Lagrangian £ {A)dS written in 
terms of (A)dS intrinsic fields is in one-to-one correspondence with the ambient-space one 
£ Amb • More precisely, considering a single term in the Lagrangian, the two descriptions 
are related by 

£ Amb ($,d$,c><9$,...) = (f) A £ (A)dS (99,V^Wv?,...), (5.1) 

where A is a constant depending on the spins and the /i-values of the fields as well as on 
the number of derivatives entering £ Amb . Regarding the action, the first attempt would be 
to consider 

/ d d+1 XC Amh =( I" dR (f ) d+A ) x ( I d d xV^C WdS ) . (5.2) 

J \J0 / \J(A)dS / 

However, the latter contains a diverging radial integral so that controlling its gauge invari- 
ance becomes ambiguous. A way of solving this problem would be to introduce a cut-off 
in order to regulate the radial integral, or similarly, a boundary for the ambient space. 
Then, the presence of the boundary breaks gauge invariance which can be restored only 
by adding boundary (total-derivative) terms in the action. The latter are the analogue of 
the Gibbons-Hawking- York boundary term needed in order to amend the Einstein-Hilbert 
action in manifolds with boundary. 

Another equivalent way is suggested by the fact that the ambient space can be consid- 
ered as a tool to rewrite intrinsic d-dimensional (A)dS quantities in a manifestly 5(9(1, d)- 
covariant form. In this respect, with the aid of a delta function, one can simply rewrite 
the (A)dS action in the ambient-space language as 

S = J Af^£ (A)dS = J d d+1 X5(V7x^-L) £ Amb . (5.3) 

As a candidate for the Lagrangian £ Amb , one may think to use the flat d-dimensional one 
where all the d-dimensional quantities are replaced by (d+ l)-dimensional ones. However, 
in general this way does not lead to a consistent (A)dS action. The reason is that, because 
of the delta function, total-derivative terms in £ Amb no longer vanish but contribute as 

6 (V7X^ -L)d xAI (■■■) = - 5' (V7X^-L)^M= (...)^0. (5.4) 

Thereby, in order to compensate these terms, the Lagrangian £ Amb has to be amended 
by additional total-derivative contributions. It is worth noticing that the latter vertices 
contain a lower number of derivatives compared to the initial vertices in £ Amb . Actually, 
this is the ambient-space analogue of what happens in the intrinsic formulation: the re- 
placement of ordinary derivatives by covariant ones requires the inclusion of additional 
lower- (covariant) derivative vertices in the Lagrangian. 

As previously discussed, a consistent (A)dS action consists of vertices containing terms 
with different number of derivatives sized by proper powers of L~ 2 : 
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In the ambient-space formalism, it is convenient to rather size such contributions by dif- 
ferent derivatives of the delta function: 



5M(V7xt-L) 5^(R-L) = {±J l j n S(R-L) , (5.6) 



since the latter naturally appear in the terms (5.4) that need to be compensated. Indeed, 
thanks to the following identity: 

«W (R -L)R*= { ~J L ^! 2]n S(R-L)R\ (5.7) 

arbitrary powers of L~ 2 can be always absorbed into derivatives of the delta function. 
Therefore, the ambient-space Lagrangian can be expanded as 

5(y^G - L) £ Amb (L" 2 ) = Y, ^ (v^* - L) C l :L , (5.8) 

n>0 

where the £ Alb's do not involve any power of L~ 2 . In order to conveniently handle the 
above series, it is useful to express 5 M by means of an auxiliary variable 5 as 



&*(R-L)=e*p(% S(R-L) [e^) g=Q . (5.9) 



For simplicity, in the following we work with the rule 8 M (R — L) = 5(R — L) (e5/L) n . 9 
The advantage of introducing the auxiliary variable 5 lies on the simple rule in dealing 
with total derivatives: 

5(V7x2 - l) d x u (•••) = - s(V7x^ -l]{x m {■■■). (5.10) 



Moreover, it also allows one to factorize the delta function in the series (5.8) and rewrite 
the Lagrangian as a polynomial function in 5/L : 

^Amb = -£-Amb(x) • (5-11) 

6 Construction of HS cubic interactions in (A)dS 

In this section we present the solution to the Noether procedure at the cubic level of for 
arbitrary symmetric HS fields in (A)dS. The logic is the same as in the flat-space case 
discussed in Section 3, thus in the following we mainly focus on those points wherein the 
peculiarities of (A)dS arise. 

Apart from the presence of the delta function, the discussions which led to the most 
general form of the TT parts of the cubic vertices still hold. The only subtleties are related 
to the total-derivative terms in (3.5) and (3.6) which no longer vanish. However, as we 



9 The 1/L in the definition of 8 has been introduced to provide a well-defined flat limit: the corresponding 
rule in flat space becomes 8(N ■ X) 8 x m {■■■) = — S(N ■ X) S N M (•••)• 
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shall explain below, their contributions can be reabsorbed into redefinitions of the cubic 
vertices. Hence, the most general expression for the TT parts of the cubic vertices reads 



[S (3 V = I d d+1 X 6(V7xl - L)cr iA2A3 (hYi,Zi) x 
x $ A i(X u Ut) $ A z(X 2 , U 2 ) <S> A *(X 3 , U 3 ) 



x i= x ' 



(6.1) 



where the Y^s and the Z^s are defined analogously to (3.13). Let us remind the reader 
that, as we have discussed in the previous section, the inhomogeneity of the vertices in the 
number of derivatives is encoded in the 5/L-dependence of the function CJ^^ . 

Whenever a gauge field joins the interactions, the cubic vertices are constrained to 
satisfy the gauge compatibility condition (3.12) associated to that field. Assuming the first 
field to be (partially-)massless (i.e. fit € {0, 1, . . . , si — 1}), one gets 

CZ^ii^z) , (U 1 -d Xi r +1 ] \ Ui=Q *0, (6.2) 

where means again equivalence modulo the d\. 's , djji • dx^s and djj.'s , i.e. modulo 
the ambient-space Fierz system. Aside from the higher-derivative nature of the gauge 
transformations, the key difference with respect to the flat case is the non-triviality of the 
total-derivative terms arising from the commutations of U\ - dx 1 with the 1^'s and the Z^s. 
Let us sketch how these total-derivative terms can be dealt with: 

• Because of the identity (5.10), the total-derivatives terms give rise to contributions 
of order 5/L and proportional to the operators X ■ &Xi or X ■ d\j i . 

• Appearing right after the delta function, the latter can be replaced by Xi ■ dxi and 
Xi ■ diji respectively. 

• Pushing these operators to the right and making them act on the fields, one can use 
the HT conditions (4.1) to replace X{ - dxi by the corresponding homogeneity degrees 
and Xi ■ djji by zero. 

All in all, one can recast the condition (6.2) into a higher-order partial differential equation 
of the form: 



ft [Yidz, - Y 3 d Z2 + | (y 2 cV 2 -Y 3 dy 3 - ^- aB ) d Yl 



n=0 



cZ A2 JhYZ)=o. 



(6.3) 

It is worth noticing that the masses of the other two fields, \ii and [i 3 , enter the equation as 
effective masses, [i^ — IYi dy 2 and n 3 — 2Y 3 dy 3 , dressed by number operators. Therefore, 
even in the massless case (//2 = = 0) a mass-like term survives. Again, depending on the 
number of (partially-)massless fields involved in the interactions, one can have a system of 
(up to three) differential equations given by the cyclic permutations of eq. (6.3). 
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7 Solutions of HS cubic interactions in (A)dS 



In this section we discuss the polynomial solutions of the system of PDEs given by eq. (6.3) 
and possible cyclic permutations thereof. Indeed, since the generating function &(X, U) is 
a formal series in U M , the latter are the only relevant ones. Our discussion mainly focuses 
on the interactions involving three massless fields which are of capital importance due to 
their connections to VE. 

7.1 Three massless case 

In the three massless case (/ij = 0, i = 1, 2, 3), one has a system of three second order PDEs 
of the form: 



Y i+1 d Zi _, - F,_! d Zi+1 + | {Y i+1 dy i+1 - Yi_i dy^) d Yi \ CZ 2a3 (i;Y,Z)=Q, (7.1) 

where [i ~ % + 3] . The latter can be solved via standard techniques (the Laplace transform 
and the method of characteristics), and its solutions are given by 

CZ 2 a 3 = exp { - 1 [Z X dy 2 dy 3 + Zx Z 2 dy 3 d G + cyclic + Z x Z 2 Z z d%] } 

xJC aiaaa3 (Yx,Y 2 ,Y 3 ,G), (7.2) 

where fC ai a 2 a 3 is an arbitrary polynomial function of the Y^s and G = Y\ Z\ + Y 2 Z 2 + Y 3 Z 3 . 
Notice that in the flat limit, one recovers the coupling (3.17). The exponential function 
provides the correct lower-derivative tails needed for the consistency of the corresponding 
(A)dS interactions. For instance, considering the lowest-derivative 4 — 4 — 4 interaction, 
JC = At^ G A , one gets 



5 . 



A 



4 



G 4 - 12 I Z x Z 2 Z 3 G 2 + 12 (|) 2 Z\ Z\ Z\ 
G 4 + ±^±£) Zx Z 2 Z 3 G 2 + 12 ( d+ L 4 ^ +6 > Z 2 Z 2 Z 2 



(7.3) 



where, in the second line we have used the identity (5.7) in order to replace the powers of 
S/L by those of L -2 . 

It is worth mentioning another way of presenting the solution (7.2). It consists in 
encoding all the 5 contributions into total derivatives as 

C IIa 2 a 3 = K-a ia2 a 3 {Yx, Y 2 , Y 2 ,G) , (7.4) 

where the Y^s and G are the (A)dS deformations: 

3 

Yi = Yi + on d ur 8 x , g = ^(y, + ft du t • d x ) Zi , 

i=l 

(ax,a 2 ,a 3 ; &) = (a, -2±1 ; ^ _ |±i, _ 2=2) , (7.5) 

of the flat-space building blocks Y^s and G . The equivalence between the two represen- 
tations (7.2) and (7.4) of the cubic interactions can be shown carrying out the integration 
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by parts of all the total-derivative terms present in (7.4). Notice that the freedom of a 
and f3 reflects a redundancy in expressing the building blocks in terms of total derivatives. 
Finally, let us conclude the discussion on the interactions of massless HS fields providing 
the example (7.3) in terms of ambient-space tensor contractions: 



[5 (3) ] TT = Xt 4 ^ j d d+1 X 5(V7X2 - L 



<Wq d M d N d p d Q $ RSTV $ RSTV + 8 § MNPQ d M d N d p $ RSTV d v 3> RSTQ 
+ 6 $> MNTV d M d N $> PQRS d R d s $ PQTV + 12 d s $ MNR T d M d N $> PQRT d R $ PQRS 



(7.6) 



7.2 General cases 



The interactions of three massless fields represent a subclass of the interactions one can 
envisage depending on the values of the /Vs • Let us notice that in the general cases the 
solutions are given by intersections of the solution spaces of the PDE (6.3) and its cyclic 
permutations. Therefore, we start our discussion from the solutions of one PDE, for which 
it is instructive to first analyze the corresponding equation in flat space (3.14). The latter 
exhibits a singular point in correspondence of the value = m3 ■ Indeed, aside from this 
value, a rescaling of m| — m| is tantamount to a rescaling of y\ . Therefore, any polynomial 
solution with 7772 7^ ms can be smoothly deformed to a solution with m-i = 777,3 , while the 
opposite is not true. Consequently, the solution space with 7772 = m% is always bigger than 
(or equal to) the one with 777-2 = m 3 ■ Indeed, an explicit analysis shows that the 7772 7^ 7773 
solutions IC(y2, y3, h-2, 113, z±) can be always expressed in terms of the 7772 = 7773 solutions 
fc-{yi-, 2/2, 2/3, 9, z\), while the opposite is not true. 

The latter phenomenon has a richer counterpart in (A)dS, where the role of tttJ — 7773 
in eq. (3.14) is played by the combinations: 

III + fi 2 - A*3 - 2 (Y 2 dy 2 -Y 3 dy 3 + n) [n = 0, . . . , m] , (7.7) 

in eq. (6.3). Indeed, because of the number operator Y2 dy 2 — Ys8y 3 , eq. (7.7) may have 
several vanishing points for \i\ + [12 — ^3 £ 2 Z . More precisely, in correspondence of the 
latter values, one can consider an ansatz of the form: 



L/ a 1 A 2 A 3 \L^ Y ^ Zj - r 2,3 A 2 -4 3 \ L ' ' 



(7.8) 



where we use Y2 for M > and Y3 for M < . Plugging this ansatz into the original 
equation (6.3), one ends up with an analogous equation for C^ A ^ A , whose re-th factor 
coincides with the operator appearing in the massless case (7.1). Therefore, the solutions 
of the massless equation provide solutions of the original equation through the ansatz (7.8). 
Notice that, when \x\ + fi2 — /J 3 ^ 2Z , the aforementioned solutions are no longer available 
since they become non-polynomial. In all the cases which can not be covered by the ansatz 
(7.8), the solutions can be expressed as arbitrary functions of the building blocks: 

Hi = • d Xl+1 d Ul+1 • - d Xl+1 ■ d Xl ^ Zi , (7.9) 
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which are the (A)dS deformations of the flat-space building blocks hi (3.16). It is worth 
noticing that this pattern is similar to what happens in flat space where the /ij-type so- 
lutions exist independently on the mass values, while the massless-type ones (involving g) 
only appear for particular values of the mi's. 

Moving to the cases in which more than one equation is involved, one has to consider 
intersections of the corresponding solution spaces. Since in flat space the only enhance- 
ment point arises for m; = m^+i , one is led to five different cases: (1) three massless 
(mi = r?i2 = = 0), (2) two massless and one massive (mi = ni2 = , 777,3 7^ 0), (3) one 
massless and two massive with different masses (mi = 0,m2 7^ "73), (4) one massless and 
two massive with equal masses (777,1 = , 777,2 = m 'i), (5) three massive. On the other hand, 
due to the presence of a richer pattern of enhancement points (/ij + /ij+i — /ij_i G 2Z), 
more combinations appear in (A)dS. The analysis of the above cases goes beyond the 
scope of the present letter, and we refer to the forthcoming paper [40] for the detailed 
discussion. 
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